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Abstract 

A Newtonian matrix cosmology, corresponding to the Banks, Fis- 
chler, Shenker and Susskind model of eleven-dimensional M-theory in 
the infinite momentum frame as a supersymmctric (0+1) M(atrix) 
model is constructed. Interesting new results are obtained, such as 
the existence of (much sought for in the past) static solutions. The 
possible interpretation of the off-diagonal entries as a background ge- 
ometry is also briefly discussed. 
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One of the most interesting recent results in string theory Q is that the 
strong coupling limit of type IIA strings is some unknown eleven dimensional 
theory, whose low energy limit is N=l supergravity in 11 dimensions. 

Although this result, combined with the general framework of string 
dualities, allows renewed hopes on the dream of unification, the fact that the 
eleven dimensional M-Theory does not have any known moduli associated 
to a coupling constant makes further progress difficult. 

One of the few attempts to describe in a concrete way M- Theory is the 
M(atrix) model of Banks, Fischler, Shenker and Susskind [§, H). The start- 
ing point is the seminal idea of || on the general description of a system of 
a given number, N, of D 0-branes by means of the dimensional reduction of 
super- Yang-Mills in (1+9) dimensions to (1+0) dimensions, i.e. a supersym- 
metric quantum mechanical model with bosonic degrees of freedom given by 
nine N by N (hermitian) matrices in the Lie algebra of U(N), X 1 . To be 
specific, the bosonic part of the action is: 

S = j dt H^D^f + ^X^X^f (1) 

where R c is the proper size of the spacelike circle (becoming lightlike when 
e | 0; this is a technical point discussed by Bigatti and Susskind in Q) and 
l\i is the scale set by the eleven-dimensional Newton's constant, G\\ = l\ v 
The bold hypothesis of BFSS is that precisely this lagrangian (when su- 
persymmetrized) describes the 11-dimensional dynamics of M-theory, albeit 
in a frame boosted in the x 11 direction, the infinite momentum frame, IMF. 
(widely used to describe some aspects of the physics of partons in the past 
|U). It has been proven that it reproduces low energy supergraviton scat- 
tering ||, and also membranes, following earlier suggestions in ||. In order 
to be able to consider arbitrary sets of particles, it seems necessary to take 
the limit N —* oo. 

It is not clear what precisely the physical meaning of the fact that D 
0-branes are described by "coordinates-matrices" is. On the one hand, it is 
clear that the eigenvalues of the matrices (wich will be denoted throughout 
by y l a , where a = 1, . . . , N ) label the ordinary positions, whereas the off- 
diagonal terms (denoted generically by z l ab ) are related to the interactions 
between the particles themselves^]. It could be that Nature is described 

1 This position seems to be at variance with the one expressed in |^|, in which the 
relative distance between two separate "clusters" in a block-diagonal matrix, with blocks 
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at some fundamental level by non-conmutative geometry, and that non- 
separability is built in the theory at a fundamental level. 

At any rate, it seems clear that the model captures at least some relevant 
aspects of interesting eleven dimensional physics. 

2 

In the IMF relativistic invariance is hidden, and only galilean one is manifest. 

k 2 

Energies are given up to a constant by H = , where the "Newtonian" mass 
is related to the eleventh component of the momentum in the IMF, ku = 
-^e"'', where uj is the hyperbolic rotation angle determining the boost (uj — > 

1 k 2 

oo). The energy in the unboosted frame is clearly given by po = ^('7+ 2^)) 
and the unboosted eleventh component of the momentum is p\\ = ~ 

k 2 

2^-). Another point of view, pioneered by Susskind, and briefly alluded to in 
the prevoius section, starts from discrete light cone quantization, in which 
the light-like coordinate x~ is compactified to a circle of circumference 2ttR. 
This means that the spectrum of p_(= rf) is discretized: p- = N/R, and we 
assume that N (to be identified later on with the number of partons) is a 
non-negative integer. Eventually the limit N — > oo is necessary. 

This means that if we want to study cosmology in this frame (the purpose 
of the present paper), it would be some form of Newtonian cosmology. It 
could be thought premature to speculate in this direction, but in our opinion 
it is always useful to be able to imagine a cosmological scenario compatible 
with the best candidate to date for a fundamental theoryQ. 

Newtonian cosmology is a beautiful theory. It is almost a deductive 
one; starting from the cosmological principle (postulating that the universe 
presents the same aspect from every point except for local irregularities) and 
assuming further that there is a universal force (the cosmological constant) 
proportional to the distance, all results of the theory flow smoothly [Q|. 

As has been stressed by Bondi, the main subject is the study of the 
motion of the substratum, which for our purposes will be idealized as the 
streaming of a uniform fluid. The concrete implementation of the cosmolog- 
ical principle will stem from the assumption that to all comoving observers 

Ya (N a x N a ) and Y b (N b x N b ), with N = N a + N b is R a , b = |*^ - but it is 

nevertheless the simplest way of getting gauge invariant expressions for the coordinates 
themselves, a very desirable property from our point of view. 

2 Some previous attempts in the broad framework of Kahiza-Klein cosmology are to be 
found in @. 
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the model presents the same appearance (at a given Newtonian time). 

Although this is the standard formulation of the cosmological principle, 
given the fact that the object of interest in the present work is a Newto- 
nian theory with extra dimensions, some modifications can be thought of. 
It could be applied to the ordinary three spatial dimensions only, for ex- 
ample, which would constitute a sort of "broken symmetry" phase. From 
this point of view, the present work is restricted to the symmetric phase. 
(Non-symmetric phases would be related to the study of cosmological (time- 
dependent) solutions of ordinary super symmetric Yang-Mills field theory (in 
n+1 dimensions), at least for toroidal T n compactifications, according to the 
standard procedure ||). Another window this sort of considerations would 
open up is the cosmological study of the effects of T-duality [14| in this 
context. 

The cosmological principle implies striking consequences on the allowed 
velocity field for the substratum. Let us consider two fundamental (that 
is, comoving with the substratum) observers, O and O', measuring the ve- 
locity of a given point P. Observer O will report V l (X J ), and observer 
O' will report V' l (X' J ). Uniqueness of the velocity at P (plus the Newto- 
nian law of addition of velocities) then necessarily implies that V H (X' 3 ) = 
V l {X 3 ) — V l (X :1 — X' 3 ) Now, the cosmological principle imposes V'(x) = 
V(x), which means that the velocity field must be linear in the coordinate la- 
bels: V t {X 3 ) = AX l B, where A and B are N x N, possibly, time-dependent 
matrices. 

We are assuming here that the substratum, is constituted somehow of 
D-0 branes, in a regime in which a matrix description of them is compul- 
sory. A fascinating question, which has been left out for future work, is the 
construction of a dual, complementary description in terms of membranes. 

In the standard scalar case, isotropy further requires that there is a 
single function of time, f(t) such that v = f(t)f. If we want this property 
to be preserved on the eigenvalues of the coordinate matrix X, and further 
demand hermiticity, we are left with 

V i (X J ) = f(t)SX i S + (2) 

where S is a unitary matrix. In the scalar case, this equation can be further 
integrated to r = R(t)fo, with / = R/R. A similar Ansatz in the present 
situation would be 

X 1 = R{t)UX l U + (3) 
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again, with a unitary U, yielding 

X* = ^X i + RU[U- 1 U,X t ]U- 1 (4) 
it 

Assuming that the stress-energy from all matter of the Universe is dy- 
namically negligible^], the only remaining force is the cosmological constant 
F l = XX 1 , (introduced in the newtonian context by Neumann and Seeliger 
in 1896 fll|| ) which means that the dynamical Newtonian equations following 
from the BFSS action principle (after a trivial rescaling) are 

x { + Yy x *> \ x ^ xi \\ = xxi ( 5 ) 

3 

The preceding equations convey an interesting generalization of the stan- 
dard vacuum equations of scalar Newtonian cosmology. The eigenvalues of 
the coordinate-matrices obey exactly the standard equations as in (Nev- 
ertheless, if the unitary matrix U is non-trivial the labels of the eigenvalues 
suffer a time-dependent permutation). The most conservative interpretation 
of the other entries of the matrices is that they are related to the interactions. 
More on this later. 

Although it is not the purpose of the present work to be exhaustive about 
it, there are very interesting new solutions. One of the simplest stems from 
the assumption that the unitary matrix U = 1, leading to 

RXi + R s Y,l x L [X j ,X*]] = ^ x o (6) 

3 

The trace of the above is 

(R - XR)trX l = (7) 

There are then two possibilities: either trX^ ^ 0, in which case the standard 
scalar result is recovered: R = Rq cosh(\/A(t — to)) (if A > 0), or R = 
Ro cos(V— A(t — to)) (if A < ); or else trX^ = 0, which means that Xq 
can be expanded in generators of su(n): Xq = ^2 a Ui a T a . To be specific, 
assuming the simplest case N = 2, and Ui a = 9(3 — i)uj a ,with furthermore 
J2 c u ic u jc = 2<5ij (equivalent to choosing a set of 3 orthonormal ordinary 
three- vectors ),then the Ui a are subject to no further restriction, whereas R 
obeys: 

R - \R + R 3 = (8) 

3 This assumption could easily be bypassed 
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whose solutions are, generically, given in terms of elliptic funcions. 

There is a fascinating particular solution, though, in the case A > 0, 
and it is a static one, R = y/\. This is very curious^, because (scalar) 
Newtonian Cosmology was abandoned essentially since the first attempts in 
the nineteenth century because it does not allow static solutions, and this 
was exactly everybody's prejudice until Milne and MacCrea revitalized the 
subject (in 1934) drawing on insights from relativistic cosmology as well as 
on Hubble's data on the expansion of the Universe. 

Apart from the static solution, R = constant, we can multiply the above 
equation with R and find 

R 2 - XR 2 + ^R 4 = C , (9) 

where C is some constant. One one can follow [11] in his analysis of solutions, 
but for the sake of briefness let's give the explicit expression for the case 
A = 0: Integration imposes that C > 0, and we find that the solution is 
given by JH 

R(t) = asd[V2a(t-t ) -) , (10) 



where a 4 = C. 

The above solution can be easily generalized to arbitrary N. Note that 
the solution is based on the asumption that we have matrices Xq such that 

^[xiAxlxiw = xl (ii) 



If we now decompose the X's on a base of su(N), using the Lie-data 
{T a , f a b c }, we see that the above equation reads 

X^J^X^X^f^fedc = X\, (12) 



where we have used the Killing metric on su(N), g a b say, to lower the index. 
It is then easy to see that iff we take the X lb such that, reverting to vector 
notation, 

J\#(#) T = 0, (13) 

3 



4 Although perhaps not so surprising, given the fact that our building blocks for the 
substratum are D-0 branes, which are well known to be BPS states, and, as such, able to 
survive in static configurations of equilibrium 
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we automatically satisfy the above equation, due to the definition of the 
Killing metric. Note that we are using Hermitean generators, so that the 
Killing metric is positive definit. 

A completely different solution is obtained for one-dimensional motion 
(still with N=2, for simplicity). A rotation can always be made so that 
Xo = 03 . Making the general ansatz U = u° + iua, it follows easily that 

/ 1-2 sin 2 9(t) - 2 sin 9 (t) cos 9 (t) 

\ -2sin0(t)cos0(t)e-^*) 2sin 2 0(t)-l 

The simplest case would be 9 = 7r/4, giving 

R - R<p 2 - XR = (14) 

Rip + 2R<p = (15) 

Again, a curious thing about these equations is that they allow static solu- 
tions when X < 0; (j) = y/—Xt + 4>o and R = Ro arbitrary. 

Given the fact that the matrix U is not trivial anymore, there is now a 
time-dependent permutation of the labels of the two fundamental particles 
in our oversimplified model. 
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The problem of the physical interpretation of the off-diagonal terms must 
now be tackled. The sort of solutions presented here (at least when A = 0) 
correspond to particles as free as they can be in the model. The most 
conservative approach would be to identify a spacetime metric such that 
their geodesies coincide with the trajectories of the eigenvalues in M(atrix) 
Cosmology. There are problems, however. The hamiltonian of the M(atrix) 
model yields an expression for the energy, which, when the effects on the 
eigenvalues y l a are separated from the rest (z l ab ) is , in a somewhat symbolic 
notation, 

e = \ jz^af + g2£4 + v(vl 4c) (16) 

a=l 

where the g l !L are quadratic functions of the y l a . This means that the classical 
equations of motion constitute a coupled system of N 2 differential equations 
of the type: 
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where the functions z % ab (t) obey: 

z: b = G* ab (yl,yl,y k c ,z l de ) (18) 

In the limit N — > oo we can assume the we have an oo 9 of fundamental 
observers, a = a 1 £ R 9 , such that for each point there is a unique trajectory 
passing through it. There are many curved spaces such that the trajectories 
Va = Haft) are geodesies of it; it suffices to take the congruence of fundamen- 
tal observers as timelike lines; parametrized by x % = a*, (so that the change 
of coordinates is defined by x l = a l (y k ,t)) and a metric, for example such 
as 

ds 2 = dt 2 - R{t) 2 5 ij dx i dx j (19) 

(which reduces to ds 2 = dt 2 on the trajectories by using ^ki/ k + Q? = 0). 
In the initial coordinates it reads: 

ds 2 = (1 - R 2 a 2 )dt 2 - R% 3 ^dy%l - 2R%^d/dt (20) 

There is a certain latitude as to how to choose the function R(t). The 
simplest way would be to demand the we recover newtonian cosmology on 
the constant time hypersurfaces (that is, on the 9-space orthogonal to the 
trajectories of our fundamental observers). This leads to its identification 
with the same function R(t) of preceding paragraphs. 

Neither should constitute a big surprise the fact that the above metric 
is not flat in general. D-0-branes only feel the metric and the RR-vector, 
so that the background constructed as above need, a priori, not satisfy the 
IIA-Sugra equations. A legitimate question is whether, building on the data 
obtained from M(atrix), we can find a IIA-Sugra solution a posteriori. 

Another point is that this technique gives naturally a 10-dimensional 
background. It could be thought more natural to interpret it as an 11- 
dimensional one, but it does not appear easy in our framework (other than 
definining 11-dimensional decompactifications of the 10-dimensional mani- 
fold M). 

4 

Quantum effects in the IMF are usually simpler to compute, because vacuum 
effects are absent. This is due to the fact that as the boost parameter uj — > oo 
those processes corresponding to diagrams which have internal lines carrying 
negative rj will go to zero, and only those with positive r\ survive. Quantum 
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cosmology would then be particularly transparent in this frame, in which 
the physical interpretation of the wave function of the universe could follow 
the guidelines of hadronic wavefunction in parton physics. We hope to be 
able to return to this topic in the future. 

It is perhaps worth stressing, to conclude, that the set of ideas here in- 
troduced with a motivation based on M-Theory, could also be contemplated 
as a natural generalization of Newtonian Cosmology per se. The basic equa- 
tions could easily have been derived just after the introduction of the matrix 
concept by Sylvester around 1850. 
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